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Program No. M-105b 
13 February 1961 
Programmer: M. V/artel 



Program Title: Modified Bessel Function of the First Kind I (x), Floating Point 
(Relocatable) 

1. PURPOSE 

1.1 To compute in floating point the modified Bessel. function I. (x). 

2. RESTRICTIONS 

2.1 Require AN-044 exponetial subroutine. 

2.2 Argument must be a positive number. 

3. METHOD 

depending 
3.1 Truncated power series or expc&ertial approximation, /on size of argument. 



3.1.1 x<7 1 {< x ^% + kC i C ~ ^ 



x„ 2k + 1 

k 



3.1.2 x£7 I «^i CI- L lL k 

*<•"* TC= 1 (Cx) 

3.2 Accuracy within *o.01% 

3.3 Reference 

H. B. Dwight, Tables of Integrals and Other Mathematical Data, 
MacMillan Co., 1947, Equations 813.2 and 814.2 

4. USAGE 

4.1 The argument must be in the A and R registers as a nomalized floating point 
number. The result is left in the A and R registers upon exiting. 

4.2 Calling sequence 

4.2.1 If AN-044 is located at (I + 150.0) 

TRA I 

Return 

4*2.2 If AN044 is elsewhere in memory 

SIX 



TRA (I + 140.0) (octal) 
P ZE EXP 
**• " f:1 - Return 
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-2- 



4.3 



4.4 



4.5 

4. 

5.1 



EXF is the first location of AN-044 

I is first location of I routine. 

After you have used 4.2.2 on a run the program can now be 
entered by 4.2.1 

Use' location 0000 to 0144 relocatable 

Temporary storage 0134 to£>137 

Error Codes 

If the limits of AN-044 are exceeded there will be a halt and transfer 
to (I, + 0044.1). 

If you press start you will return to the main program. 

Both L and V loops are used by this routine. 

Coding Information 



Constants (floating point) 



1 

2TT 



M. 



M 



*3 



0100 
0102 
0104 
0106 
0110 
0112 
0114 
0116 
0120 
0122 
0124 
0126 
0130 
0132 



7.0 
1.0 
2.0 

6.283185308 

12.0 

144.0 

2880.0 

86400,0 

3628800,0 

203212800.0 

3.0 

7.5 

52.5 

8.0 



XlOOCU0--tJ-5^QO2£.^Q___S-AX Jhlosl 

+0C0075.O ADD 13 February 1961 



0001 .0 +42.0044.1 STA" 
_-f^Q0pp_._0___XA_R. 



jQQ02^.Q_4:&)_c0077 t Q___STO. 
+30.0075.0 FCA 



0003.0 +06.0100.0 FSB 
+_52»QQk5.D___TBP_ 



-00a^0_+A.QaL0 J .O.__CTL. 
+66.0020.0 CTV 



0005.0 +30.0076.0 FCA 

4--a^03Xi4.a__EDV_ 



.^QQ6« r Q-+-3-^7-7?2*a.--ESJ- 
+07.7772.0 FMP 



0007.0 +35.77T6.0 FST 

+-57.r7-76o.a--TBA_. 



001 . +Q 7*77 7 g *Q FMP 



+35.777^.0 fst 



~66i"i .0 +05.61 o"4.o Tdv 
+P.4 8 .77.7-2 i Q__F_Ap. 



0012.0+^.7772.0 FST 
+30.7776.0 FCA 



00 T 3". "+07 ."7774" :o" " FMP" 
____„___±3^777ji.?___F_ST 



0014.0 +05.0110.0 FDV 
+04.7772.0 FAD 



0015.0 +35.7772.0 FST 
±_\QaiZ£ t 0__F_CA_ 



0016.0 +07.777JU0 FMP 
+35.7774.0 FST 



OOI7.O +05.0112.0 FDV 
_-_______+a*U77-72.0__Ji:AD_ 



-0020^Q-+35-777-^0-P-SI. 
+30.7776.0 FCA 



0021.0 +64.00^0.0 CTL 
+57.77h0.0 TRA 



-0022^0 _+ocu.oooa«.o_ _ -_-_- 

-00.0000.0 — 



0023.0 +00.0000.0 — - 

_ „ „ _ •qo..qqqg*q- --.v. 



.4Q2K&-+QMQ0QJQ— «™ Pamr* 

-oo.oooo.o — M ^ 05b 

002570 TbbV666b;6""~--- IS" February- 196l 

- O O.QQQO tQ — 



QQZS^Q „+QO^QOOO ._0 _ „ _~ - - . 
-00.0000.0 

0027.0 4-06.666b.~6~~--- 
-00.0000.0 — 



-0030.0- 4--Q? .777*U& - -EMP. 
4-35.777^.0 FST 



0031 .0 +05.011^.0 FDV 

+0S»7-772*a__EAD.. 

-0032.0-. 4--35 .7-772 . Q. _ .EST. 
4-30.7776.0 FCA 



0033.0 +07.777^.0 FMP 

.003k .0- *o§ . au.6. a _ _edy_ 

4-04.7772.0 FAD 



OO35.O 4-35.7772.0 FST 
^30. 7X76. Q_ . _EC A. 



.oa36»o-.d-n7.zr7^*a __emp_ 

4-35.777U.0 FST 



0037.0 +6H.oo^o.o cTT 
_„.±5.7.7l60*D— TBA. 



OOUO-0 4-0S.0120.0 FDV 



4-04.7772.0 FAD 



00*11.0 4-35.7772.0 FST 

4-„3Q«77-76,0._FCA. 



00^2.0 4-07. 777*4.0 FMP 



4-35.777^.0 FST 



00^3.0 4-05.0122.0 FDV 

.___: — 4--04,7772^a-FAa. 



OQ l OUQ +35.7770 . FST 



)7. 



4-57.OOOO.O TRA 



OOH5.O 4-30.0076.0 FCA 
-4-57.030O.0___TRA. 



00U6.0 4-00.777^.0 — CIA 



4-77.00^.1 HTR 

00V7.0 +6h.b&p.b ctl 
4-^7.77£Q-0___tea. 



-QQ*)O-.-e-+6&fQ96&r0— G-TV Page~5~ 

+35-7776.0 FST M . 105 b 



0651 76 +Vd."6ib6.6" fca 

+07.0076.0 FMP 
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jX52_jQ _*_^777k*.Q- - FJ-T. 
+447777^.0 FSQ 



0653.O "+V5". /777~4 •" 0" " FSf " 
+30.7776.0 FCA 



_005k-0.=hQ5.7-77^-Jl— EDV. 
+35.OI3H.O FST 

OO55.O -R0.0102.~0 FCA 
4-05,0132.0 — EX1V_ 



-0056, 0- +05 , 0074»O- -FDV- 
+35.777^.0 FST 



0057.0 +07.012*1.0 FMP 

13547776,0 — £S£- 



00^0,0- 4-~30^7-77^0l__FCAc. 
+07.777^.0 FMP 



0061 .0 +35.0136.0 FST 
+0.7*QJ2_6 S .0__J : :MP_ 



_Q0£a.0_±pi t 77_l6 ft 0__.FAp. 
+35.7776.0 FST 



OO63.O +64.0070.0 C~T 
+5_7-776O i 0__TRA 



oo64 % o_+p_o.qppo. P___--~ 
-00. 6000. — 



0065.0 +00.0000. r* — 

r_0Q_OPOj3.O___-.--_ 



._Q0M_Q_+_0_Q_QQp_Q_0___---_ 
-00.0000.0 



0067.0 +00.0000.0 

_-l3Q.0QQQ.Q___--- 



0070. Q 4.30.0136.0 — EXA 



+07.777^.0 FMP 



0071.0 +07.013Q.0 FMP 

__ 4-a^777-6--0— ""-AD. 



0072 . +35 , 7776 . FST 



+30.0102.0 FCA 



OO73.O +06.7776.0 

.________4-a7--oi-3-^0-- 



FSB 

JJ_P. 



■ovjk.-o -+5T-. ooW- J TRA _. . 

+00.0000.0 CLA j£g^ 

OO^Vo"+^r6o^\"6" ---""" 13"»teuiiy"19eL" 



,00-^0000^1 ---. 



0076.0 +00.0000.0 



-00.0000.0 



0077.0 +00.0000.0 — 

. -- - =,00-00.00.0. - - - . 



-04 00.0- j-_70«0QQ0Jl— ==■_- 
-00.0000.0 — 

0101.0 +00.0000.0 

- »QQ. 0001 .1 



O402.Q.-+ikQ.QOQ0*O.__ARS.. 
-00.0000.0 — 



0103.0 +00.0000.0 — 
-00.0000.1 m 



-Oiak»0-+_ta*QQQQ..JCL_-ARSL 
-00.0000.0 — 



0105.0 +00.0000.0 

-00.0001 .0 — 



oj-o6*o._+_62,. 2077-0— ---- 

+52.42H0.0 TRP 



0107.0 +00.0000.0 — 
-00.0001 .1 =. 



.OUO.OL +60,OQQ0^0__._SJ0 
-00.0000.0 — 



01 11 .0 +00.0000.0 — 

_-Q0 J .QQQ2-Jl-.--rr_ 



-04J-2.0_+M.OQQ0^O.__ESQ.. 
-00.0000.0 — 



oi 13.0 +00.0000.0 

— -____ J oo,OQaU,.a-_™.- 



.04 4-^.0- +-55*0000^0- — tsc 
-00.0000.0 — 



0115.0 +00.0000.0 

00^0006.0----— . 



Q14-6.-0- .+52.liOO*0— IBP. 

-00.0000.0 



0117.0 +00.0000.0 — - 
^00.001Qa_------ 



01 2a.G +67* 2-7&)«0 CFV 

-00.0000.0 Page 7 
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=QQvQ013-JD--~- 



0122.0 +fr). 3^62.0 STO 



-00.0000.0 - 



0123.0 +00.0000.0 

rrQOJXMjLQ— ----. 



0124.0 +fr). 0000.0 SICL 



■00.0000.0 — 



0125.0 +00.0000.0 

^oq.oqou-Q— ==•_-_ 



0126,0 4-7?UQQQQ.Q — BNC 



-00.0000.0 



0127.0 +00.0000.0 

G0.QGQU4— «- *- 



-00.0000.0 



0131 .0 +00.0000.0 

-00-0003.0 



01 ia^O_+iiO t OOOO.O___ARS 

~~-oo.6ooo.~o --- 



0133.0 +00.0000.0 

-00.0002.0 



QL3^._o +P_q._opoo.p___— - 

-00.0000.0 --- 



0135.0 +00.0000.0 

-00.0000.0 — 



0136 t o.+op. poop. 0___—- 
-00.0606.6 — 



0137.0 +00.0000.0 

-00.0000.0 =ri 



.-oito,.o-+4-5.oa76*o— sax. 

+^2.01^.1 STA 



01 41 .0 +01 .01 Ml. ADD 
+115^00^^.1— 5TA- 



_0 1 +2-*Q -+to.O00Q^0— ARS - 
+00.0000.0 CLA 



OIU3.O +42.0045.1 STA 
+5.7,0001-^1— -TRA. 

-ov4iM>-+oa^oooa^o— -_— . 

•.nn.nnm .n 



